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Abstract 



It is shown that the nuclear shadowing of charmonium due to the modification 
of the nuclear parton distribution is similar in the factorization approach based 
on non relativistic QCD and in the color evaporation model. In the first model, 
a separate study of the color octet and color singlet contributions to the yields of 
the various charmonium states as well as the contributions of these states to the 
total J/^S production is performed. It is found a clear xp dependence of these 
contributions which can reproduce experimental data for moderate xp. 
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1 INTRODUCTION 



Since the observation of a possible anomalous suppression in lead on lead collisions by 
the NA50 collaboration [I], ^ a lot of theoretical work has been done trying to explain the 
data. Some papers use conventional physics ||-|| but most of the explanations are based 
on QGP formation [3|-[]l0|, the new state that has been predicted to be formed at high 
enough temperatures and densities of hadronic matter and from which J/\I/ suppression 
has been proposed as a signature long ago fill. Although there are some consensus in 
the possibility of this state to be reached in 158 AGeV/c lead on lead collisions, a lot of 
uncertainties remain to be solved on the production and suppression of the J/\l/. A deep 
study of the absorption mechanisms seems necessary in order to assert the anomaly on 
the observed suppression. 

Charmonium production in hadron-hadron collisions is now quite accurately described 
by a factorization approach (FA)[] based on non-relativistic QCD (NRQCD) fL2f , al- 
though the so called color evaporation model (CEM) has also a remarkable success in 
describing integrated as well as differential cross sections [13| [ I4[ . The main difference 



between them is the treatment of color, being the CEM quite radical in this sense: it 
simply forgets about color. In the calculations of this paper both models will be used, 
turning out to give almost identical results. 

In parton model, the production mechanism is described by the fusion of a parton 
coming from the projectile with longitudinal momentum fraction x% and a parton from 
the target with longitudinal momentum fraction x 2 - For a fixed xp > ( X p = X\ — ^2/: 
this momentum fraction are X\ ~ Xp and X2 — M 2 jx\S. 

For collisions involving nuclei the situation is much more complex. Multiple scattering 
of the initial partons as well as of the produced pre-resonant states with nucleons from 
the nucleus and even of the produced resonances with other produced particles (co- 
movers) must be taken into account. It is very common to call initial state interactions 
to multiple interactions of the initial partons (mostly gluons in the case of J/\I/ and 
and final state interactions to the rescattering of the pre-resonant state. These names 
are senseless at high energies, where time ordering is lost. It is possible, however, to 
factorize the multiple scattering of light partons from that of the heavy pre-resonance. 
The former can be described by modifying the nucleon structure function inside nuclei 
(EMC effect). For the latter a probabilistic formula for the heavy system to interact is 
usually employed. This formula can also be derived in a field theoretical approach. All 
these features are discussed in |f[5f| . In this way the situation is different in open charm 
as well as in Drell-Yan (DY) production, where there is no absorption of the produced 
heavy system (at least at moderate energy) but the effect of nuclear modifications on 
structure functions is almost the same as in the case of charmonium [I5|. A combined 



1 Sometimes it is called Color Octet Model, although it seems now more frequent FA. The name COM 
is not very correct because it is not a model, but the solution of an effective theory and not only color 
octet states are taken into account. 
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study of the different reactions should help to distinguish the contributions of each effect 
to the suppression. 

In this paper we are only interested in studying the effect of the modification of the 
nucleon structure functions inside nuclei (EMC effect). As we will see, this modification 
has a big influence on the absorption of charmonium at large xp, a problem that remains 



unsolved |L7j . Although this effect is negligible in the integrated cross section at present 
energies, it will be important at higher ones. A better understanding of the different 
suppression mechanisms in pA and AB collisions could shed some light on what is hap- 
pening in PbPb, and, for that, the Xp dependence of absorption will be very helpful. 
Concerning suppression, large xp is related with higher energies in two senses: first, 
from the point of view of suppression by EMC effect, as the momentum fraction of the 
target parton that takes part in the collision is x 2 — M 2 /sxf, larger xp as well as s 
imply smaller x 2 (the region of shadowing corrections). Second, from the point of view 
of rescattering of produced pre- resonances by nuclear matter (not studied here), again 
larger ^ or s means faster produced particles, that could lead to a bigger absorption 
cross section [I8[ This is discussed in detail in reference flT5| |. 



Further uncertainties come from the treatment of color in absorption of pre-resonant 
state by nuclear matter. It is usually assumed that a pre-resonant color octet state is 
absorbed by rescattering with other nucleons inside nuclei. But, as it is well known by 
FA, J/ty is not fully formed in color octet state and, as we will see, the fraction of the 
different color contributions to the production cross section is affected by the EMC effect, 
being the color singlet states less shadowed than the color octet ones. Then, although 
the correction by EMC effect itself is not very big, it will have some effect in the amount 
of absorption by nuclear rescattering if different absorption cross section are used to color 
octet and color singlet states as done, for example, in reference [|J. 

2 QUARKONIUM PRODUCTION CROSS SECTION 

The production cross section for a quarkonium state H according to the process 

P + T -> H + X 
can be described in framework of parton model by the integral 

a H = dx 1 dx 2 ^2f l P (x 1 )ff(x 2 )cr(i,j-*H), (I) 
J i,j 

where f i \X\) is the i-parton distribution function in the projectil (target). 
The FA based on NRQCD gives the partonic elementary cross section, 



^3-*H) = J2C% Q[n] <0»>, (2) 
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in terms of short distance parts Cqqua describing the partonic process to obtain a QQ 

pair in a state [n] and long distance coefficients < > which describe the hadronization 
of this state into a physical particle H. The short distance parts can be calculated in 
series of a s {2mo). However, the long distance ones are non-perturbative quantities, so 
they are not exactly calculable at present. They can be taken as parameters to be fixed 
by experimental data or lattice calculations. The matrix elements for the production of 
charmonium and bottonium states were computed in [p0fl -|p2|. 

One of the quantum numbers which describe the state [n] is the color of the cc pair. 
Then, in this approach, contributions coming from different color states of the pre- 
resonant cc pair are separately taken into account. 



An alternative model for quarkonia production is the CEM [|13|. 14]. In this model 
it is assumed that the different quarkonia states come from the QQ pairs created below 
the threshold for open charm production (i.e. s = x\x 2 s = Am 2 D ). In this case the cross 
section for the production of a given particle H is assumed to be a constant fraction Ffj 
of the hidden cross section Gqq defined as: (U 



a QQ\ 

so that 



4m D 



ds \ dx\dx 2 f i (xi)fj (x 2 )a(s)8(s — xix 2 s), (3) 

4m? 



v p n-+h(s) = F p H d QQ {s), (4) 



being Fjj the parameter for each particle H . This parameter is free and independent of 
the process. a(s) in equation (|) is the partonic cross section to produce a pair cc with 
and invariant mass s calculated in pQCD p3 |. 



For J/ty production, contributions coming from the different decays must be taken 
into account. At fixed target energies the main contributions are the decays of the 
various Xcj states and that of the \E r/ . Other contributions, like the decay of B states are 



negligible ]20[]. Thus, 



<Tj,9 = off 9 + B(*'^J/*)<T< lt ,+ B (XcJ - J/*)<T XeJ , (5) 

7=0,1,2 

being crfL, the direct J/^f production, i.e. not coming from decays, and B(H — ► J/^f) 
the branching ratios for particle H to decay into a J/ty. Experimentally it is found that 
the direct J/ty contribution is about 60% of the total J/ty cross section, decay is less 
than 10% and all Xcj contribute with more than 30%. 

Decays into ^' are not important, and then ty' production is dominated by direct 
production. 

In CEM it is not necessary to make this distinction, because the decays are implicitly 
taken into account in the parameters Ffj. The FA on the contrary offers the possibility 
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to study these different contributions in a separate way, calculations in [^(|-[^2j give an 
accurate description of these ratios. 

Another important point is the contribution of the different color states to the pro- 
duction of these particles, i.e. the color content of the pre-resonant state. In fact, the 
FA gives that direct J/^f production is almost completely produced in color octet state, 

is also predominantly (about 90% or more) produced in color octet, and the main 
contribution of the Xcj states to J are given by color singlet matrix elements. Then, 
a separate study of the EMC effect in different particles and color states is possible in 
this approach. As we will see, the color octet and color singlet contributions to the 
production of charmonium have different suppression by EMC effect. 



To introduce nuclear corrections by the EMC effect, we use the results of In this 
reference a parametrization of experimental data for the ratio 

Aff(x, Q 2 ) 

is evolved by means of DGLAP equations. For this parametrization only DIS and DY 
experimental data has been taken into account. In this formula, A is the atomic num- 
ber, ff^(x,Q 2 ) is the i-parton distribution function for nucleons inside nucleus A and 
ff{x,Q 2 ) is the corresponding one in a free nucleon. Thus, the structure functions of 
equation ([!]) must be substituted by the corresponding ones for nucleons inside nuclei: 
A/f(x,Q 2 )i?f(*,Q 2 ). 

These ratios have a shadowing region for x < 0.1, in which the parton distribution 
function inside nucleus is depleted {Rf < 1) in comparison with that of free nucleon, an 
antishadowing region for 0.1 < x < 0.3 in which this ratio is bigger than 1, a second 
shadowing region for large x and a second antishadowing region (usually atributed to 
Fermi Motion of nucleons) for x > 0.85. 



3 NUMERICAL CALCULATIONS 



The cross sections in the FA have been calculated using CTEQ-LO |25[ parton distribu- 
tion functions, m c = 1.5 GeV and /i = 2rriQ, i.e. the same parameters as in |2t| in order 
to use the same values for the matrix elements < > fitted by them. In the case of 



CEM we use GRV-HO parton distribution functions and the same values for the c 



quark mass and renormalization scale. 

The results are presented in the usual parametrizationf] 



a pA(AA) = A a a pp_ ^ 



2 It is worth noting that this common parameter a is in fact dependent on A and then, although 
useful, it must be taken with caution. 
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It is clear that CEM will give the same results for all the different charmonium states, 
because the factors Ffj are canceled when the ratio is done. It is not so in the case of 
FA, where every particle has different contributions coming from the different matrix 
elements and in the case of J/\l/ also from the decays of the other charmonium states. 
In Fig.l it is shown the a parameters for the production of total J/^f, direct J/^, 
and the contribution of the various xcj states to the J/^f in proton-gold collisions. 
Different color contributions are presented in separate curves, being each curve the result 
of (^) in which the ratio is done taking the same contribution of color (i.e. Vgtigiet/ a singled 

'octet / a octet an d a total / a total) ■ We can see here the general features of the nuclear effects 
in charmonium production. First of all, the results for the CEM and for the FA are 
not so different, and difference is much smaller than the absorption produced by other 
mechanisms, like nuclear absorption of pre-resonant states by nuclear matter or co- 
movers. Second, shadowing corrections are more important at higher energies, as it is 
expected, because predominant contribution to the cross sections comes from partons in 
the target with xi ~ M 2 /s, so as energy raises, we are going to lower and lower values 
of x, the region of larger shadowing. Finally, color singlet states are less shadowed than 
color octet ones (color singlet curves are always above color octet ones). Moreover, more 
than 90 % of direct J/\& is produced in color octet, while, when different charmonium 
decays into this particle are taken into account, this ratio is less than 70%. This gives 
direct J/ty production more suppressed than total one. In Fig. 2 we present our results 
for gold on gold collisions. The same general trends as in proton-gold are present in 
this case. We can observe that effects now are comparatively less important than in the 
proton- go Id case. 

We have also calculated the xp dependence of nuclear effects for these processes. Figs. 
3 and 4 are the a's for proton-gold and gold-gold collisions for ^fs = 39GeV. We find 
that shadowing corrections are much more important as xf is getting larger, as expected 
(see above). This fact is found experimentally. Direct comparison of our results with 
experimental data is not possible because rescattering of produced pre-resonances is the 
main contribution to absorption. However, we are allowed to compare how much of the 
experimentally found bigger absorption at large Xp is due to shadowing in the structure 
functions. To do this, we consider rescattering corrections to be independent on ip. If 
both effects are taking into account in separate exponents 



a pA = A *(x P ) A (3-l a pp^ (g) 

where a(xp) are our results and (3 is a new absorption parameter independent on 
and parametrizing other effects different from EMC, we can compare our results with 
experimental data by normalizing our curves in this way to one of the less shadowed 
experimental points. This has been done in Figure 5 where experimental points are 
from reference ||27Hp|. This normalization corresponds to a abs ~ 6 mb in the probabilistic 

3 In this reference different targets were used and a parametrization in the form (m) is given taken 
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formula for nuclear absorption, which is smaller than a abs ~ 7mb, used for instance in 
and [|nj. A different shape is found for large enough xp (not reached for being 
the experimental data more absorbed. Then EMC effect alone can not be responsible 
for this suppression, the additional depletion coming from other effects like an increase 
of rescattering absorption at large xp, intrinsic charm, etc. The study of the EMC effect 



in charmonium production has been done previously in references f28| , |29| , |30|j . In this 



references, scaling in x 2 of the ratios a pA /a pp is found for some energies. This scaling 
is ruled out in reference p7| . In our approach this scaling violation should be due a 
different final state absorption depending on if and s. Reference [|^| restrict their 
study to the < xf < 0.55 region. In this sense, our results are in agreement with 
theirs. The main differences between the present work and the previous ones are the 
use of charmonium production cross section taken from hadron-hadron collisions and 
parametrization of EMC effect from DIS and DY experimental data. In particular, in 



29| the parametrization for the ratio (|6]) is much more shadowed than the one used 
by us for low x and the a abs used there to reproduce experimental data is quite small 
(a abs =2.6 mb). 



4 CONCLUSIONS 

Effects of the modification of nucleon parton distribution functions inside nucleus on 
the production of different charmonium states have been studied. Energy as well as xp 
dependence of this effects have been calculated for proton-gold and gold-gold collisions. 
We have found a depletion on the ratio between the pA(AA) to pp cross sections for 
higher energies and larger xf- This effect is, at present energies, much less important 
that the other absorption mechanisms of charmonium production, like pre-resonant ab- 
sorption by rescattering with other nucleons in the nucleus or by co-movers. However, 
these processes superimpose to the production mechanism, so they are affected by the 
corrections due to EMC effect. 

When FA is used, a separate study of the contributions of different color states and 
particle decays to the production mechanism is possible. Our result is that color octet 
states are more shadowed than color singlet ones, and that a slight difference is obtained 
for J/ty production when decays of other charmonium states are taken into account. 

We have also found a slightly different absorption for J/\l/ and states, in the case of 
FA. This difference is not present in experimental data, being too small to be seen over 
the other absorption mechanisms. We can say then that nuclear corrections for J/\l/ and 

production due to EMC effect are almost the same. In the CEM, as we have already 
said, they are exactly the same. We also obtain that CEM and FA give almost the same 

into account all the different nuclei. As we have said, this parameter is in fact A-dependent, then for 
comparison with our calculations only data for W target is used (we supose that nuclear corrections to 
PDF in Au and in W are almost the same). 
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results in all our calculations. 

Finally, we obtain a different shape for experimental and theoretical curves in if, 
being the suppression in experimental data bigger than that in our calculations for large 
xf- Although our curves almost fit the experimental points (within errors) for the xf 
dependence of nuclear suppression, data seems to prefer a bigger absorption for large xp. 
So we can conclude that modification of the nuclear parton distribution functions can 
account for the xp dependence of J/^ suppression for moderate xp and that other mech- 
anisms apart from EMC effect is needed at large xp in order to reproduce experimental 
data. 
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Figure captions 



Fig.l. Center of mass energy dependence of a for p-Au collisions. Different lines are: 
FA total contribution (solid), singlet contribution (dotted), octet contribution (dashed) 
and CEM (dotted-dashed) 



Fig. 2. Same as Fig.l but for Au-Au collisions. 



Fig. 3 xp dependence of a for y/s = 39 GeV p-Au collisions with the same conventions 
for lines as in Fig.l. 



Fig. 4. Same as Fig.3 but for Au-Au collisions. 



Fig. 5. Comparison of xf dependence of nuclear corrections by EMC effect with ex- 
perimental data of 800 GeV/c protons incident on tungsten target from reference [27] 



Theoretical calculations have been normalized to one of the less shadowed experimental 
points (see text). Lines follow the same convention as Fig.l. 
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